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ABSTRACT 

Orthogonal  Frequency  Division  Multiple  Access  (OFDMA)  is  proposed  in  IEEE  802. 16e  standards  for 
broadband,  high  data  rate  fixed/mobile  communications.  As  the  basic  modulation  technique  in  OFDMA  is  orthogonal 
frequency  division  multiplexing  (OFDM),  it  is  severely  affected  by  carrier  frequency  offset  (CFO).  The  CFO  causes  loss  of 
orthogonality  among  the  subcarriers.  In  this  paper,  we  propose  a  method  for  estimation  of  CFO.  The  proposed  method  is 
based  on  the  preamble  suggested  for  OFDMA  mode  of  WMAN  where  the  preamble  is  loaded  on  every  third  subcarrier. 
The  proposed  method  is  an  extension  of  Schimdl  and  Cox's  method.  Simulations  are  conducted  for  AWGN  and  frequency 
selective  channels  to  demonstrate  the  performance  of  the  proposed  method.  CRLB  for  the  proposed  fractional  CFO  is 
derived  and  is  shown  that  the  variance  of  proposed  CFO  closely  follows  CRLB. 
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1.  INTRODUCTION 

An  OFDM  system  effectively  consists  of  N  sinusoidal  sub  carriers  with  frequency  spacing  l/T,  where  T  is  the 
active  symbol  period  of  each  subcarrier.  The  k'h  subcarrier  will  thus  be  at  fk=f0+klT,  where  f0  is  a  reference  frequency. 
Without  loss  of  generality,  we  can  assume  /o=0.  The  modulated  subcarriers  overlap  spectrally,  but  since  they  are 
orthogonal  over  a  symbol  duration,  they  can  be  easily  recovered  as  long  as  the  channel  does  not  destroy  the  orthogonality. 
An  unwindowed  OFDM  system  has  rectangular  symbol  shapes.  Hence,  in  the  frequency  domain  the  individual  subchannels 
will  have  the  form  of  sine  functions.  The  sine  function  is  defined  as 


sinc(x) = 


Kx 


(1) 


With  this  definition,  it  can  be  confirmed  that  sinc(0)  =  1  and  that  zero  crossings  occur  at  +1,  ±2,  +3,  ....  Since  the 
sine  waves  existing  in  each  OFDM  symbol  are  truncated  every  T  seconds,  the  width  of  the  main  lobe  of  the  subcarrier  sine 
functions  is  2/T,  i.e.,  there  are  zero  crossings  every  1/rHz.  and  the  first  side  lobe  is  only  13  dB  below  the  main  lobe  of  the 
subcarrier  [1,2]. 

Therefore,  N  subcarriers  can  be  packed  into  a  bandwidth  of  NIT  Hz,  with  the  tails  of  the  subcarriers  trailing  off  on 
either  side.  Since  the  zero  crossings  of  the  frequency  domain  sine  pulses  all  line  up,  as  long  as  the  frequency  offset  =  0, 
there  is  no  interference  between  the  subcarriers.  One  intuitive  interpretation  for  this  is  that  since  the  FFT  is  essentially  a 
frequency-sampling  operation,  if  the  frequency  offset  is  negligible,  the  receiver  simply  samples  received  signal  at  the  peak 
points  of  the  sine  functions,  where  the  ICI  is  zero  from  all  the  neighboring  subcarriers. 

In  practice,  of  course,  the  frequency  offset  is  not  always  zero.  The  major  causes  for  this  are  mismatched 
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oscillators  at  the  transmitter  and  the  receiver  and  Doppler  frequency  shifts  owing  to  mobility  [3].  Since  precise  crystal 
oscillators  are  expensive,  tolerating  some  degree  of  frequency  offset  is  essential  in  a  consumer  OFDM  system  such  as 
WiMAX.  For  example,  if  an  oscillator  is  accurate  to  0.1  parts  per  million  (ppm),  /offset25  (/c)(0-l/>pm).  If  3.6  GHz  and  the 
Doppler  is  fd=400Hz,  /offset-  (/c)  (Q.lppm)  +  fd  =  760Hz,  which  will  degrade  the  orthogonality  of  the  received  signal, 
since  now  the  received  samples  of  the  FFT  will  contain  interference  from  the  adjacent  subcarriers  [4]. 

Research  has  shown  [5]  that  the  frequency  offset  error  must  less  than  4%  of  the  OFDM  subcarrier  interval  in 
AWGN  channel  and  less  than  2%  of  the  OFDM  subcarrier  interval  in  Rayleigh  channel  to  satisfy  the  requirement  of 
system  synchronization. 

If  the  CFO  is  several  times  larger  than  the  subchannel  spacing,  the  offset  can  be  split  into  an  integer  part 
(an  integer  multiple  of  the  subcarrier  spacing)  and  a  fractional  part  i.e.,  "e  +  p  where  p  represents  the  integer  part  of 

normalized  CFO  and  e  represents  the  fractional  part  of  normalized  CFO.  The  integer  part  of  the  normalized  CFO 
introduces  a  shift  in  the  subcarrier  position  at  which  true  data  is  available.  The  fractional  part  of  CFO  attenuates  and  phase 
rotates  the  symbol  on  the  subcarrier,  and  also  causes  inter  carrier  interference  [5].  All  these  effects  of  CFO  leads  to  very 
high  bit  error  rate  (BER),  sometimes  even  half  of  the  total  data  symbols  may  be  in  error  [6],  To  summarize,  OFDM  is 
severely  affected  by  carrier  frequency  offsets.  Thus,  to  mitigate  the  effects  of  frequency  offset  between  the  transmitter  and 
the  receiver,  we  need  to  have  a  very  efficient  frequency  synchronization  technique. 

CFO  estimation  techniques  may  be  classified  as  time-domain  (pre-FFT)  or  frequency  domain  (post-FFT) 
techniques.  Time-domain  methods  are  generally  used  to  estimate  the  fractional  part  of  the  CFO,  although  some  of  these 
techniques  can  also  estimate  the  integer  part.  The  post-FFT  techniques  are  usually  used  to  estimate  the  integer  part  of  the 
CFO  (normalized  by  the  subcarrier  spacing)  after  the  fractional  part  has  been  identified  and  corrected. 

1.1  Existing  CFO  Estimators 

In  order  to  mitigate  frequency  offset  effect,  various  techniques  have  been  proposed  to  estimate  the  CFO  for 
OFDM  systems.  The  frequency  offset  mitigation  using  null  symbols  is  proposed  in  [7]  but  the  use  of  a  null  symbol  is  not 
practical  in  a  burst  mode  system  as  the  transmission  will  always  be  a  null  while  no  data  is  transmitted. 
The  self-cancellation  schemes  were  presented  in  [8]  and  [9],  where  the  same  data  is  sent  in  more  than  one  subchannel. 
These  schemes  are  not  very  practical  as  they  reduce  the  number  of  useful  subchannels  by  at  least  half.  In  [5], 
Moose  proposed  a  maximum  likelihood  (ML)  estimator  using  repeated  data  symbol.  Although  the  symbol  values  need  not 
be  known,  repetition  of  OFDM  symbols  is  in  essence  tantamount  to  a  training-based  scheme  as  it  utilizes  extra  bandwidth. 
Data-assisted  frequency  acquisition  and  tracking  were  proposed  in  [10],  where  periodically  inserted  known  symbols  were 
explicitly  used.  In  [11],  Schmidl  and  Cox  proposed  a  training  symbol-based  frequency  synchronization  that  utilized  an 
OFDM  symbol  with  identical  halves. 

This  was  later  generalized  to  a  training  symbol  with  multiple  identical  parts  [12].  Various  blind  techniques  have 
also  been  proposed.  In  [13],  Van  de  Beek  et  al.  developed  an  ML  estimator  by  exploiting  the  redundancy  in  the  cyclic 
prefix  (CP).  This  method,  however,  is  developed  based  on  the  nondispersive  channel  model  and  suffers  the  error  floor 
effect  in  the  presence  of  frequency-selective  fading  channels.  Schmidl  and  Cox  proposed  in  [14]  a  blind  estimation  method 
that  is  only  suitable  to  recover  CFO  values  that  are  multiples  of  the  carrier  spacing.  In  [15-16],  Liu  and  Tureli  took 
advantage  of  the  presence  of  virtual  carriers  in  OFDM  signaling  and  proposed  blind  estimation  methods  reminiscent  of 


Impact  Factor  (JCC):  4.9467 


Index  Copernicus  Value  (ICV):  3.0 


Carrier  Frequency  Offset  Estimator  for  OFDMA  Mode  of  WMAN 


95 


spectral  analysis  techniques  in  array  processing,  i.e.,  MUSIC  and  ESPRIT.  It  was  later  shown  that  the  proposed  MUSIC 
algorithm  is  indeed  the  ML  estimate  of  the  CFO  with  a  virtual  carrier  present  signal  model  [17,18].  The  drawback  is  that  it 
usually  requires  multiple  OFDM  symbols  to  achieve  desirable  performance,  thus  introducing  an  extra  delay  at  the  receiver. 
Tufvesson  proposed  a  frequency  offset  estimation  approach  using  repeated  pseudonoise  (PN)  sequences  [19].  Its  estimation 
range  is  large.  The  performance  of  the  time  synchronization  method  using  PN  sequence  is  very  nice.  But  there  is  conflict  in 
choosing  the  period  of  PN  sequence.  Large  period  of  PN  sequence  means  better  time  synchronization  performance,  but  this 
will  result  in  smaller  range  of  frequency  offset  estimation.  Now,  let  us  check  how  far  the  previous  techniques  are  useful  in 
frame  synchronization  of  OFDMA  mode  of  WMAN. 

1.2  Suitability  of  Existing  CFO  Estimators  in  Application  to  OFDMA  Mode  of  WMAN 

To  implement  CP  based  frame  synchronization  algorithm,  CP  length  must  be  known  to  the  receiver. 
IEEE  802.16e-2005  [20]  standards  specify  4  different  CPs  for  OFDMA  mode  of  WMAN.  They  are  64,  128,  256,  and  512. 
BS  selects  one  CP  and  transmits  the  data.  As  CP  is  not  known  to  SS,  CP  based  frame  synchronization  algorithm  has  to  be 
repeated  for  all  four  possible  CPs,  which  is  not  only  time  consuming  but  also  laborious  task.  Also,  these  algorithms  are 
suitable  only  for  AWGN  and  very  slow  fading  channels.  But,  BWA  channels  based  on  NLOS  propagation  are 
characterized  by  the  stringent  multi  path  delay  profile  and  the  fading  characteristics.  Hence,  CP  based  frame 
synchronization  algorithms  are  not  suitable  for  OFDMA  mode  of  WMAN. 

Preamble  based  frame  detection  algorithms  depend  on  the  structure  of  preamble  used.  Most  of  these  algorithms 
are  based  on  preamble  having  two  (or  four)  identical  parts,  these  can  be  generated  by  loading  on  every  second  (or  fourth) 
subcarrier  of  the  preamble.  But  IEEE  802.16  standards  [20]  specify  loading  on  every  third  subcarrier  of  the  preamble  with 
FFT  size  2048  (or  4096)  i.e.,  getting  three  identical  parts  is  not  possible.  Hence,  preamble  based  frame  detection  algorithms 
too  are  not  suitable  for  OFDMA  mode  of  WMAN. 

Therefore,  we  have  to  develop  a  frequency  synchronization  algorithm  especially  for  OFDMA  mode  of  WMAN. 

This  paper  is  organized  as  follows.  Section  2  presents  the  proposed  method  for  estimation  fractional  CFO. 
CRLB  for  the  proposed  fractional  CFO  is  derived  in  Section  3.  Section  4  discusses  the  estimation  of  integer  frequency 
offset.  Simulation  results  are  provided  in  Section  5. 

2.  PROPOSED  METHOD  FOR  FRACTIONAL  CFO  ESTIMATION 

A  PN  sequence  is  transmitted  along  with  the  OFDMA  preamble.  For  each  segment,  there  are  32  preamble  time 
series  [20].  96  PN  sequences  are  specified  for  these  preamble  time  series.  One  out  of  these  96  known  PN  sequences  is 
transmitted  by  BS  and  the  SS  has  to  identify  the  PN  sequence  transmitted  and  cell  ID  used.  Cell  ID  represents  a  unique 
number  corresponding  to  the  PN  sequence  used  by  the  BS.  Maximum  likelihood  estimation  [21]  can  be  used  to  find  the  PN 
sequence  transmitted.  The  procedure  is  as  follows:  At  the  receiver,  N  samples  are  collected  from  frame  boundary  of  the 
received  signal  sequence.  These  N  received  signal  samples  are  correlated  with  each  of  the  96  known  PN  sequences. 
When  the  transmitted  PN  sequence  is  aligned  with  receiver  PN  sequence,  a  correlation  peak  is  observed.  In  a  multipath 
channel,  the  largest  correlation  peak  occurs  at  the  energy  peak  of  the  multipath  delay  profile.  Hence  frame  boundary 
aligned  to  the  strongest  multipath  gives  the  largest  correlation  peak.  The  main  drawback  with  this  method  is  that  it  fails 
with  the  presence  of  frequency  offset.  It  should  be  applied  after  frequency  compensation.  The  uncompensated  phase 
difference  between  signal  samples,  as  a  result  of  frequency  offset,  will  destroy  the  correlation  property  of  the  PN 
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sequence.  We  have  developed  a  simple  algorithm  [22]  which  can  be  used  for  joint  estimation  of  PN  sequence  and  fractional 
frequency  offset.  In  this  algorithm,  the  PN  sequence  correlation  is  robust  to  frequency  offset  in  the  absence  of  noise. 

Let  the  time  domain  samples  corresponding  to  the  preamble  and  data  symbols  be  denoted  by  xp(ri)  and  xd(n), 
respectively.  Then,  the  transmitted  signal  x(ri)  can  be  written  as 

x(n)  =  x  (n);     0<n<N  +  L-1 

(2) 

=  xd{n)\     n>N  +  L 

The  internal  block  correlation  properties  of  the  preamble,  discussed  in  following  subsection,  are  used  in 
development  of  this  algorithm  to  estimate  fractional  frequency  offset  and  cell  ID.  The  preamble  carrier-sets  are  defined 
using  the  following  formula  [20] 

k  =  n  +  3.m  (3) 

Where  k  denotes  indices  of  subcarriers  allocated  to  a  specific  sector,  n  is  the  index  of  the  carrier-set  number 
(ranging  from  0  to  2)  and  m  is  a  running  index  from  0  to  567. 

Therefore,  each  sector  eventually  modulates  each  third  subcarrier.  Figure  1  shows  three  carrier  sets  corresponding 
to  three  sectors.  For  easy  identification,  different  colors  are  used  for  different  carrier  sets.  Subcarrier  index  0  corresponds  to 
the  first  used  subcarrier  of  the  preamble  (Figure  1).  A  total  of  568  subcarriers  at  the  indices  0,  3,  6,  9,. ..,1701  form  the 
carrier  set  for  sector  0.  Similarly,  a  total  of  568  subcarriers  at  the  indices  1,  4,  7,  1702  form  the  carrier  set  for  sector  1 
where  as  a  total  of  567  subcarriers  at  the  indices  2,  5,  8,  1700  form  the  carrier  set  for  sector2.  The  used  subcarriers  in  the 
preamble  are  modulated  using  a  BPSK  modulation  with  a  specific  PN  sequence.  The  selection  of  PN  sequence  depends  on 
the  sector  used  in  the  cell  [20].  The  subsymbols  for  the  guard  carriers  and  dc  carrier  are  zero  amplitude  symbols. 
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Figure  2:  Time  Domain  Structure  for  the  Preamble  in  OFDMA  Mode  of  WMAN 

Figure  2  shows  the  block  structure  for  maximum  correlation  of  496.  Let  us  denote  three  blocks  of  682  samples  Au 
Bh  and  C\  representing  a^n),  b\(ri)  and  c^n)  respectively  as  follows. 

a1(n)  =  xp(n  +  L);       n  =  0,1,.  ..,681 


bl(n)  =  xp(n  +  L  +  6S3);       n  =  0,1,...,  681  (5) 

c{{n)  =  x;)(«  +  L  +  1365);      n  =  0,1,.  ..,681  (6) 

2.1  Proposed  CFO  Estimator 

By  following  the  procedure  similar  to  [21],  we  estimate  the  fractional  frequency  offset  and  cell  ID.  In  [21], 
Schmidl  and  Cox  exploited  the  correlation  existing  in  the  two  identical  parts  from  the  frame  boundary  of  the  preamble  for 
the  estimation  of  the  fractional  frequency  offset.  Correlation  between  two  identical  parts  results  in  energy  of  one  part. 
But,  in  the  preamble  corresponding  to  OFDMA  mode  of  WMAN,  two  identical  parts  do  not  exist.  Instead,  this  preamble 
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consists  of  three  correlated  blocks  as  given  in  (3)  to  (5).  The  correlation  between  any  two  correlated  blocks  results  in  some 
residual  phase.  So,  Schmidl  and  Cox  method  is  modified  [22],  as  shown  below,  to  suit  the  preamble  corresponding  to 
OFDMA  mode  of  WMAN. 

2.2  Preamble  Structure  Recommended  for  OFDMA  Mode  of  WMAN 

OFDM  modulation  with  2048  subcarriers  is  used  in  OFDMA  mode  of  WMAN.  The  first  OFDM  symbol  in 
downlink  of  WMAN  is  the  preamble.  The  preamble  is  generated  with  1702  used  carriers,  346  null  carriers  -  173  on  the  left, 
172  on  the  right  as  guard  bands  and  a  dc  carrier,  as  shown  in  Figure  1. 


173 

851  subcarriers 

851  subcarriers 

172 

Zeros 

BPSK  symbols  on 

BPSK  symbols  on 

Zeros 

every  3rd  used 

every  3rd  used 

carrier 

carrier 

Left 

guard 

band 


1703  used^carriers,  Nusai 
Null  symbol 


Right 
guard 
band 


Figure  1:  Preamble  Subcarriers  are  Divided  into  Null  Carriers,  Dc  Carrier  and  Used  Carriers 

The  cell  area  corresponding  to  downlink  is  divided  into  three  segments  and  these  segments  are  named  as  sector  0, 
sector  1,  and  sector  2.  The  used  subcarriers  are  divided  into  three  sets  corresponding  to  the  three  sectors  of  the  cell. 

By  multiplying  a  sample  in  the  block  B{  with  corresponding  conjugate  of  a  sample  in  block  A,,  we  form  a 
reference  sequence  used  for  finding  correlation.  From  (3)  and  (4),  the  sequence  v(n)  is  formed  as 


v{n)  =  a"l(n).b1(n);      n  =  0,1,...,  681 


(6) 


There  are  96  such  sequences  corresponding  to  96  PN  sequences  specified  for  OFDMA  [20].  Now,  the  t 
correlation  sequence,  where  ;'  =  1,  2,...,  96,  is  given  by 


v('\n)  =  a^inyb^in);     n  =  0,1,. ..,681; 


(7) 


Where  a\  (ri)  and  b\  (ri)  represent  the  samples  in  block  Ai  and  block  B[  corresponding  to  the  f  PN  sequence  of 
the  preamble.  Substituting  (3)  and  (4)  in  (7),  the  f  correlation  sequence  becomes 


v(i\n)  =  x('r(n  +  L).x('}(n  +  L  +  683);  n  =  0,1, ...,681; 


(8) 


To  keep  the  exposition  simple,  let  us  assume  an  ideal  channel  with  no  noise.  Then,  samples  of  the  received  signal 
are  given  by 


r{n)  =  elN  jx(n) 


(9) 


Here  e  represents  the  CFO  between  transmitter  oscillator  and  receiver  oscillator  normalized  by  subcarrier  spacing, 
given  by  £=f/Af  where  in  Af  is  the  subcarrier  spacing  and  /  is  the  carrier  frequency  offset,  80  is  an  initial  arbitrary  carrier 
phase,  and  u(n)  is  a  zero  mean  circularly  symmetric  complex  white  Gaussian  noise  with  variance  a2  ■ 
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From  the  estimated  frame  boundary  dopt,  we  find  the  sequence 

r\n)  =  r\dvt  +  n).r(dopt+ L  +  6S3);  n  =  0,1,. ..,681;  (10) 
We  then  find  the  correlation  of  all  sequences  V(l)(n),  i  =  1,2,      96  with  r  (n).  Substituting  (9)  in  (10),  we  get 

r\n)  =  x"{dopt  +n).x{dopi  +n  +  683)e  N  ;  n  =  0,1,..., 681;  At  dopt=  L,  r(ri)  contains  the  samples  of  preamble 
as  shown  below. 

,1366;z£ 

r\n)  =  xp{L  +  n).xp{L  +  n  +  6^)e'^r;  n  =  0,1  681; 

Now,  the  PN  sequence  index  corresponding  to  transmitted  PN  sequence  is  given  by  the  value  of  i  corresponding 
to  maximum  correlation  magnitude  [22]  i.e., 

681 

Rv,r(i)  =  ^v*(n).rXny,  i  =  1,2,2...,96  (11) 

And  PNseq  =  arg  max{|^r(/)|}  (12) 

Where  PNseq  is  the  PN  sequence  number.  Then,  unique  identification  number  (cell  ID)  corresponding  to  this  PN 
sequence  number  is  obtained.  Out  of  96  PN  sequences  [20],  first  32  PN  sequences  correspond  to  Segment  0,  next  32  PN 
sequences  for  Segment  1  and  the  remaining  32  PN  sequences  correspond  to  Segment  2.  Now  depending  on  the  PN 
sequence  number,  we  can  easily  identify  the  Segment  selected  by  the  Base  station. 

Next,  let  us  consider  the  estimation  of  frequency  offset.  The  correlation  between  v*(n)  and  r'(n)  with  maximum 
magnitude  is  given  by 

.1366 m  68i 

K~=e    N  X|*,(L  +  W)f  k(L  +  w  +  683)f  (13) 

Now,  the  fractional  frequency  offset  is  estimated  from  the  angle  of  ^?max  as  [22] 
N  , 

e  =  \R  (14) 

1366^"^ 

As  every  third  subcarrier  of  the  preamble  is  loaded  in  OFDMA  mode  of  WMAN,  the  maximum  detectable 
normalized  fractional  frequency  offset  becomes  \e  I  =  1.5.  For  example,  let  frequency  offset,  e  =1.7  is  added  to  the  signal  at 

the  receiver  input.  Using  this  algorithm,  we  get  the  estimated  frequency  offset,  £=—1.3.  That  means,  after  fractional 
frequency  offset  compensation,  we  observe  a  subcarrier  rotation  of  3.  This  subcarrier  rotation  is  considered  as  integer 
frequency  offset. 

The  variance  and  Cramer  Rao  Lower  Bound  (CRLB)  of  the  proposed  CFO  estimator  are  derived  in  the  next 

Section. 
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3.  STATISTICAL  ANALYSIS  OF  THE  CFO  ESTIMATOR 

Theoretical  expressions  for  Variance  and  Cramer  Rao  Lower  Bound  (CRLB)  are  derived  as  shown  below. 

3.1  Variance  of  the  CFO  Estimator 

Assuming  perfect  timing  synchronization,  we  have  estimated  the  frequency  offset  using  partial  periodic  structure 
present  in  OFDMA  preamble.  It  is  observed  from  the  expression  of  r'(«)  that  a  pair  of  samples  of  two  consecutive  periodic 
sections  after  correlation  with  transmitted  preamble  data  from  a  noise  free  channel  will  experience  a  phase  offset  of 
1366K/N.  So,  phase  offset  between  a  pair  of  periodic  sections  is  obtained  from  its  complex  conjugate  multiplication. 
Complex  conjugate  multiplication  of  nth  sample  and  (n  +  683)th  sample  can  be  expressed  as 

r  \n)  =  r'(n  +  djr{n  +  683  +  dopt )  (15) 

When  complex  conjugate  of  first  sample  of  useful  part  of  transmitted  preamble  is  multiplied  with  its  683rd  sample 

we  get 

v(n)  =  x'p(n  + L)xp(n  +  6S3  + L)  (16) 
Considering  perfect  time  synchronization  i.e.,  dopt  =  L,  the  correlation  between  r\ri)  and  v(n)  is  given  by 


(17) 


z(n)  =  V  (n)r'(n) 

=  xp(n  +  L)x'p(n  +  683  +  L)r{n  +  L)r{n  +  683  +  L) 

Let  us  consider  AWGN  channel.  Then  the  received  signal  can  be  written  as 

7.x 

r(n)  =  e ' N"' \p{n)  +  v{n)  (18) 

Where  e  is  the  normalized  frequency  offset,  and  v  is  the  nth  time  domain  AWGN  sample  added  to  OFDMA 
preamble  with  mean  zero  and  variance  ov2.  Substituting  (18)  in  (17),  we  get 


z{n)  =  \xp(n  +  L)\  \xp(n  +  L  +  683)|  e 


2    1 — ^7— 


j—e(n+L) 


+  \xp(n  +  L)\  jcp(n  +  L  +  683)y(«  +  L  +  683)e  N 

+  |x;,(«  +  683  +  L)|  xp(n  +  L)v  (n  +  L)e  N 

+x  (n  +  L)x*  (n  +  L  +  683) V  (n  +  L) v(n  +  L  +  683) 


Which  can  be  further  written  as 


1366  ire 

i2 


z(n)  =  \xp(n  +  L)\  |xp(w  +  L  +  683)|  e  N  +v(n)  (19) 
Where       is  the  effective  noise  component  of  the  complex  conjugate  product  given  by 
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  I  |2  -J — f(n+L) 

V(n)  =  \xp(n  +  L)\  xp(n  +  L  +  683)y(n+L  +  683)e  N 

+  bcp(n  +  683  +  L)|  xp(n  +  L)v  (n  +  L)e  N 

+xp  (n  +  L)x*  (n  +  L  +  683)V*(ra  +  L)v("  +  L  +  683) 

Phase  offset  estimate  "(j)(n)  from  the  nth  sample  and  (n+683)',!  sample  of  the  OFDMA  preamble  may  be  represented 

as 

*   s  1366tze 

=  — ; — +«L»  (2°) 


Where  <j>noise{ri)  is  the  phase  of  the  effective  noise  component  y(n) .  Phase  offset  <j>  estimated  from  pair  of  M 
(=  682)  samples  of  the  correlated  sections  present  in  OFDMA  preamble  is  represented  as 

~    13667K-     1  ,  ^ 

iV  Mto 

Estimated  frequency  offset  e  from  the  two  consecutive  correlated  sections  of  the  OFDMA  preamble  can  now  be 
written  as 

i  M-l 

— y  $  o) 

e  =  f.  +  iT£^!50   (22) 

1366^£-/A^ 

Noise  at  the  «th  sample  of  OFDMA  preamble  is  Gaussian  distributed  and  may  introduce  a  uniformly  distributed 
phase  offset  of  ±Jt  with  zero  mean  and  variance  7T/3.  Therefore,  variance  of  frequency  offset  from  (22)  may  be  represented 

as 

N2 

var(f)  =   (23) 

3M(1366)2 

Substituting  N  =  2048  and  M  =  682,  we  get 

var(£)  =  l.lxl0~3  (24) 

3.2  CRLB  of  Estimated  CFO 

To  derive  the  Cramer  Rao  Bound  of  the  estimated  frequency  offset,  we  rewrite  (19)  as 

Z(n)  =  a(n)  +  jj3(n)  (25) 
1366^' 


Where  a(n)  =  A  cos 


and  f}(n)  =  A  sin 


N 

1366flg' 
V     N  J 


+  V  (n)  (26) 


+  V(n)  (27) 
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Wherein  A  =  |xp(n  +  L)  \x  (n  +  L  +  683) 
Now,  we  want  to  estimate  the  parameter  vector 

r  =  [A,ef 

The  joint  probability  density  function  of  z(n)  can  be  written  as 


f{z(n);y)  =  Kexp^—2 


M-l.  , 

X \(a(n) -«,(«))  +(A«)-A(«))  ) 


where 


K  = 


1 


,  =  A  cos(0) ,  /?,(«)  =  A  sin(0) 


Where  in  0  =  1366^£-/A^  (30) 

Taking  the  first  order  partial  derivatives  of  /  (z(n);y)  with  respect  to  e  yields 


9[ln(/(z(n);r))]_  -1 
3e  2o-2 


=  2fT2  ^ 


X  ((«(n)  -  A  coS(</>)f  +  (/?(«)  -  A  sin  W)2) 
f  1366^^1 


2(a(n)-A  cos(<*))(A  sim 
+  2(a(n)-A  sin(<*))(-A  cos(<*)) 


{  N  J 
N  J 


Which  can  be  further  written  as 


a[ln(/(z(n);r))]    -1^1366*^,.,    .      ,   u  . 


(T2  TV 

+  (a(n)-A  sin((Z>))(-A  cos(^))] 

Similarly,  the  second  order  partial  derivative  of  /(z(n);y)  with  respect  to  e  yields 

32[ln(/(z(«);r))]  -lfl366^Y^V,2  .  2^,  — ,  ,^ 
—  ^  1  =  ~^  (-^)  ZLA  sin-(^)  +  vmA  cos(^) 

+  A2cos2(0)  +  vqpA  sin((Z*)J 
Which  can  be  further  simplified  as 


1(7)  = 


(28) 


(29) 


'92[ln(/(z(»);;r)]1 

'a2[ln(/(z(n);r)]]" 

H 

3A2 

H 

3A3£-  j 

>[ln(/(z(n);r)]] 

'a2[ln(/(Z(n);r)]] 

3A3£- 

1  H 

Now,  the  variance  of  the  fractional  CFO  estimator  is 
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Var(e)  >  ■ 


1(7),,  Jd2[\n{f{z{n);Y)] 


de2 


(31) 


N 


1366^-V      ,    (1366^r)2  M(SNR) 
'  MA2 


N  J 


Where  SNR  =  ■ 


(32) 


Therefore,  the  Cramer  Rao  Lower  Bound  (CRLB)  for  the  variance  of  frequency  offset  estimation  is  given  as 


(var(£))  = 

v  '  CRLB 


N1 


'crlb  (l3667ryM(SNR) 
Substituting  N  =  2048  and  M  =  682,  we  get  the  CRLB  for  the  variance  of  estimated  frequency  offset  as 


(33) 


(var(f)) 


3.3394x10" 


SNR 

4.  SIMULATION  RESULTS 


(34) 


In  our  simulation,  a  wireless  system  operating  at  3.6GHz  and  with  bandwidth  of  20MHz  is  assumed. 
Three  channel  scenarios  are  considered:  an  AWGN  channel,  time-invariant  NLOS  Rayleigh  fading  channels 
(SUI-2,  SUI-4  and  SUI-6),  and  time-variant  Rayleigh  fading  channels  (3GPP  Ped-B  and  3GPP  Veh-A)  with  the  Maximum 
Doppler  Shift  400  Hz.  The  simulation  parameters  of  OFDMA  system  given  in  Table  1  are  used.  The  preamble  is  generated 
with  subcarrier  loading  as  shown  in  Figure  1.  Every  third  used  carrier  is  loaded  with  PN  sequence  corresponding  to 
Sector  =  0  and  cell  ID  =  10  [20].  A  frequency  offset  of  1.7  times  the  subcarrier  spacing  is  assumed  in  the  simulations. 

For  this  frequency  offset,  the  estimated  fractional  frequency  should  be  e  —  —  1 .3 .  In  what  follows,  we  analyze  the 
performance  of  the  proposed  fractional  CFO  estimator. 

Table  1:  Simulation  Parameters  of  OFDMA  System 


Parameter 

Channel  Bandwidth,  BW 

20  MHz 

Oversampling  rate 
n  =  FS/BW 

28/25 

Sampling  Frequency 

Fs  =  floor(«.W/8000)x8000 

22.4  MHz 

No.  of  Total  Subcarriers, 

2048 

Subcarrier  Spacing 
A/  =  FS/JV 

10.94  KHz 

Useful  Symbol  Time 
ATS  =  l/Af 

91.43  us 

Guard  Time 
Tg  =  G.TS 

22.86,  11,43,5.7,  2.86  us 

Left  Guard  Band 

Zeros  are  loaded  on  left  173  subcarriers 

Right  Guard  Band 

Zeros  are  loaded  on  left  172  subcarriers 

Modulation  (Preamble) 

BPSK  on  every  third  subcarrier 

Modulation  (Data  Symbols) 

QPSK 
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Table  1:  Contd., 


Carrier  Frequency 

3.6  GHz 

MS  Velocity 

30  kmph  for  3GPP  Ped  B 
120  kmph  for  3GPP  Veh  A 

Max.  Doppler  Frequency 

100  Hz  for  3GPP  Ped  B 
400  Hz  for  3GPP  Veh  A 

Number  of  Cells 

3 

4.1  Performance  of  Fractional  CFO 

To  illustrate  how  the  frequency  offset  estimation  algorithm  performs  in  both  AWGN  and  frequency  selective 
channels,  we  have  performed  simulations  with  1000  different  realizations  of  channel  and  noise.  SUI  channel  models 
(for  fixed  WMAN)  [23]  and  3GPP  channel  models  (for  mobile  WMAN)  [24]  are  used  to  test  the  proposed  methods  in 
frequency  selective  channels.  fdT  =  0.04  is  considered  in  3GPP  channel  model.  The  impulse  response  of  the  channel  is 
normalized  to  unit  norm.  From  the  results  of  1000  realizations,  we  computed  mean  and  variance  of  the  fractional  frequency 
offset  estimate  considering  frame  boundary  along  the  strongest  path  of  the  channel.  The  computed  values  of  mean  and 
variance  of  the  fractional  frequency  offset  are  tabulated  in  Table  2.  Table  2  shows  that  for  AWGN  as  well  as  frequency 
selective  channels,  as  SNR  increases  from  0  dB  to  10  dB,  the  error  in  mean  value  decreases.  At  SNR  =  10  dB,  3GPP 
Pedestrian-B  channel  exhibits  largest  error  variance.  The  second  largest  error  variance  is  exhibited  by  SUI-4  channel. 
The  high  error  variance  exhibited  by  3GPP  Pedestrian-B  and  SUI-4  channels  is  due  their  bad  power  profiles  as  compared 
to  all  other  frequency  selective  channels.  Remaining  channels  -  AWGN,  SUI-2,  SUI-4  and  3GPP  Vehicular-A  channels 
exhibit  more  or  less  similar  error  variance.  Combining  the  results  of  low  error  variance  of  frequency  selective  channels 
with  the  results  of  mean,  we  observe  that  the  fractional  CFO  estimate  is  very  close  to  the  true  value. 

Table  2:  Mean  and  Variance  of  the  Fractional  CFO  Estimate  e  (the  Actual  Normalized  CFO,  e  =  1.4) 


S.  No. 

Channel 
Type 

SNR 

Fractional  CFO  Estimate 

Mean 

Variance 

1. 

AWGN 

OdB 

1.3972 

0.7189X10-3 

10  dB 

1.3988 

0.0558X10-3 

2. 

SUI-2 

OdB 

1.3957 

0.4577x1 0"3 

10  dB 

1.3985 

0.0559X10-3 

3. 

SUI-4 

OdB 

1.4027 

0.5002X10"3 

10  dB 

1.4003 

0.0629X10"3 

4. 

SUI -6 

OdB 

1.4027 

0.4714X10"3 

10  dB 

1.4002 

0.0530X10-3 

5. 

Ped-B 

OdB 

1.4066 

0.5780X10"3 

10  dB 

1.4027 

0.0952X10"3 

6. 

Ped- A 

OdB 

1.3985 

0.4644x1 0"3 

10  dB 

1.3991 

0.0507xlff3 

We  repeated  the  experiment  for  different  values  of  SNR.  1000  different  channel  and  noise  realizations  are 
conducted  at  each  SNR  value  with  SNR  ranging  from  0  dB  to  40  dB  in  steps  of  5  dB.  A  frequency  offset  of  1.4  times  the 
subcarrier  spacing  is  assumed  in  the  simulations.  Figure  3  and  Figure  4  show  the  mean  and  variance  of  fractional  CFO 
estimator  over  different  SNR  values  when  frame  boundary  is  considered  along  the  strongest  path  of  the  channel. 
These  plots  show  that  the  fractional  CFO  estimator  shows  less  deviation  of  mean  from  the  true  value,  even  in  the  presence 
of  high  doppler  channels  like  3GPP  Pedestrian-B  and  Vehicular-A  channels.  And  also,  the  proposed  fractional  CFO 
estimator  exhibits  more  or  less  similar  performance  in  all  channels  and  yields  good  results  even  at  SNR  =  0  dB. 
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Figure  3:  Mean  of  Fractional  CFO  Estimate  (Normalized  CFO,  e  =1.4,  Frame 
Boundary  is  Considered  along  Strongest  Path  of  the  Channel,  and  /Dt=0.04) 


0  5  10         15         20         25         30         35  40 

SWIindB 


Figure  4:  Variance  of  Fractional  CFO  Estimate  (Normalized  CFO,  8  =1.4,  Frame 
Boundary  is  Considered  along  Strongest  Path  of  the  Channel,  and  /Dt=0.04) 

From  Figure  4,  we  note  that  the  variance  CFO  estimate  increases  as  frequency  selectivity  increases.  The  channel 
dispersion  length  of  Pedestrian-B  channel  is  large  as  compared  to  Vehicular-A  channel  and  also,  its  second  channel  tap 
amplitude  is  largest  among  all.  This  causes  high  error  variance.  Therefore,  the  variance  of  fractional  CFO  estimator  with 
3GPP  Pedestrian-B  channel  is  more  than  that  of  3GPP  Vehicular-A  channel.  In  the  similar  way,  the  variance  of  SUT4 
channel  is  more  than  that  of  SUI-6  channel.  In  this  case  the  bad  power  profile  of  SUT4  increases  the  error  variance  of  the 
CFO  estimator.  If  we  compare  frequency  selective  channels,  SUI-2  channel,  due  to  its  line  of  sight  component,  exhibits 
least  minimum  error  variance. 


0  5  10         15         20         25         30         35  40 
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Figure  5:  Variance  of  Fractional  CFO  Estimate  in  Comparison  to  its  CRLB  (Normalized 
CFO,  e  =1.4,  Frame  Boundary  is  Considered  along  Strongest  Path  of  the  Channel) 

From  Figure  5,  we  note  that  the  error  variance  of  proposed  CFO  estimator  in  AWGN  closely  follows  the  CRLB 

plot. 
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5.  CONCLUSIONS 

In  this  paper,  we  proposed  a  method  for  joint  estimation  of  fractional  CFO  for  OFDMA  mode  of  WMAN. 
Error  Variance  and  CRLB  of  the  proposed  CFO  estimator  are  derived.  Simulations  are  conducted  for  AWGN  and 
frequency  selective  channels  to  demonstrate  the  performance  of  the  proposed  method.  It  is  found  that  the  error  variance  of 
proposed  CFO  estimator  in  AWGN  closely  follows  the  CRLB  plot.  Also,  the  estimated  fractional  CFO  is  very  close  to  the 
true  value  for  both  AWGN  and  frequency  selective  channels. 
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